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1. Introduction 



Evolution kernels are the main ingredients of the well-known evolution equations for the parton 
distribution in DIS processes [1] and for the parton distribution amplitudes [2] in hard exclusive 
reactions. These equations describe the dependence of the parton distributions on the renormal- 
ization scale //^. Previously, two-loop calculations were performed for the unpolarized forward 
DGLAP evolution kernel P{z) in [1, 3, 4] and, what is more cumbersome, for the nonforward 
ERBL kernel V{x,y) in [2, 5, 6] that was challenging and complicated technical tasks. Here, we 
present the results of a direct calculation of evolution kernels for the transversity distributions in 
next-to-leading-order (NLO) performed in the MS scheme. These calculations are carried out in 
the Feynman gauge within a single mold for both the forward kernels and the nonforward ones, 
i.e., 

P^{x) = as P^{x) +alP^{x) +... , (1.1) 
V^{x, y) = a, V^{x, y) + a] V^{x, y) + . . . , (1.2) 

where ag = as(;[i^)/(47r). 

Note that the kernel P^ was first obtained in [7] within a light-cone gauge calculation and 
shortly thereafter the corresponding anomalous dimensions 7f (n) were presented in [8, 9]. The 
kernel was reconstructed in [10] on the basis of the knowledge of the structure of symmetry- 
breaking terms for the kernel, which first appeared at the two-loop level. For the reader's con- 
venience, let us explain these issues in more detail. Those terms of that are responsible for 
the conformal-symmetry breaking can be fixed and expressed via some special convolutions of the 
known [10-12] one-loop kernel elements. At the same time, the remaining part or, in other words, 
the symmetrical part (in terms of a conformal-group representation) of this kernel can eventually 
be restored from a certain part of the forward kernel P^ [10]. This possibility of "guessing" will 
not be pursued here. 

The calculation of V'^ or P^ can be performed following the standard procedure to find the 
renormalization-group generators in the MS scheme. Expressions for them in terms of the renor- 
malization constant Zr for every diagram F in (see, e.g. , [14] and [6]) are given by 

Zr = 1 - KR'iT), V {P) = -asda, (z«) = a,a„, (i^ii?'(F)) 2k,R\T). (1.3) 

Here, (i) R' is the incomplete BPHZ i?-operation; D — A — 2e \s the space-time dimension, (ii) 
K separates out poles in £, whereas Ki picks out a simple pole, and (ill) Z^^^^ is the coefficient of 
the simple pole in the expansion of Zy- To introduce an appropriate notation for the analysis of 
the two-loop results, let us start with the leading order Vq (Pq) results obtained in a covariant 
^-gauge^ 

P^ix) = Cf [p,{xU-5{l-x)] , (1.4) 
If (x, y) = Cf [2F^(x, y)+ + {x ^ x, y ^ y)] - 5{y - x) (1.5) 

Here, Po(x) = ; F'^(x,y) = ; x = 1 — x, y = 1 — y; symbol (. . .) i denotes different 

1 — X y y — X 

distributions like p(a;)+ — p{x) — 5{l—x) p{z) dz a,ndV{x,y)+ — V{x,y) — 5{y — x) jQV{z,y) dz. 



^ The gauge parameter ^ is defined via the gluon propagator in lowest-order perturbation theory which reads 
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The diagrammatic expansion of the kernels is presented in the Table below, where ^-dependent 
terms appear in the partial diagrams a, c canceling out each other in the complete results in 
Eqs. (1.4-1.5), as expected. The slash on the line of each of these diagrams denotes the delta 

Table 1: Diagrammatic expansion of the one- loop kernels with MC denoting the mirror-conjugated 
diagrams 




function 6{x — nk/nP), where k is the momentum on this line, while n is a light-cone vector 
{n^ — 0) . These diagrammatic calculation rules can be traced to the momentum representation of 
the composite operator ■ijj{0)a^i,i/j{\n), denoted here by 0, and were elaborated in detail in [6]. The 
abbreviation MC in the figures denotes the mirror-conjugate diagrams, while the symbol C denotes 
the corresponding mirror conjugation of arguments, C9{y > x)f{x,y) = 9{y > x)f{x,y) + 9{y < 
x)f{x,y). The local current, corresponding to the operator (g), is not conserved and, therefore, 
there is no "plus" prescription imposed on Eqs. (1.4), (1.5). Therefore, the separate 5-functions 
survive. 

Noting that the product (yy) V^{x,y) is symmetric under the exchange x ^ y, one realizes 
that the corresponding anomalous dimension matrix can be diagonalized in the Gegenbauer basis 
{ipnix) = {xx)Cn'^{'2x — !)}• The deeper reason for this is that conformal symmetry survives at 
the LO level [15]. On the other hand, at NLO the conformal symmetry does not hold true (in 
the MS scheme) owing to rcnormalization effects. These generate specific terms in that break 
this X ^ y symmetry as well as the diagonalization property mentioned above. For brevity, the 
corresponding terms will be referred to as "nondiagonal (diagonal)" ones. 

In the next section, the contributions to Pi and for each of the 2-loop diagrams will be 
demonstrated explicitly. In Sec. 3, we analyze the structure of both calculated kernels which are in 
accord with the expected manifestation of these symmetry breaking terms. Finally wc confirm the 
results for P^, calculated in [7], as well as the result for found in [10]. The kernel provides 
the key ingredient, necessary for any complete NLO analysis of exclusive processes involving 
transversely polarized vector mesons via a QCD evolution of their distribution amplitudes. For an 
illustration of this NLO evolution, we analyze in Sec. 4 how it affects the transversely polarized p- 
meson distribution amplitude (DA) at the characteristic scale //^ applicable to the 5-meson scmi- 
leptonic decay [13]. Our main findings are summarized in Sec. 5 together with our conclusions. 
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2. Diagram-by-diagram presentation for and 

Here, we present the diagram- by-diagram results of the calculation of the DGLAP, , and 
ERBL, Vi ^ kernels at the two-loop level for ^ = 0. In all there are 19 diagrams in the list below 
where we also display the diagrams with a zero contribution. The full list of them can be found in 
[4] . Note that superscripts -k mark the obtained new result for each diagram. The results for the 
other diagrams can be restored from those obeying the DGLAP [4] or ERBL [6] evolution kernels. 
Diagrams f* and h with gluon-loop insertions include also the corresponding ghost loops. Let us 
remark that there are only four basic scalar topologies of integrals, the latter being presented ^ in 
[6]. 





P(x) = -C|.[po(l + lnx)]^ 


V(x,y) = -2Cl 


Ce{y > x)F^ (l + In (l - ^))]^ 




P(x) =CFTrNf^5{l-x) 


Y{^,y) =CFTrNf^5{y-x) 




16 

P(x) = -CfCa-5{1-x) 


V(x,y) = -CFCA^5(:y-x) 


g /i 




r /20 4 \i 
P(x) = -CpTrNf \po + 3 InxjJ ^ 


r / 20 4 T \ 1 

V(x,y) = -2CFT,.Nf [Ce{y > x)F^ [- + 


h /! 


'3 \ 


P(x) =CfCa[i 


/31 5, \i 

On \ — In X 

"V9 3 J\ + 


V(x,y) = 2CfCa 





^ see also corrections to these results in Appendix B in [16] 
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P(x) =-ACf(Cf-^ 



In a; 



V(x,y) = -2Cf(Cf - ^)ce{y >x)Un{l- ^) 




P(x) = ACf{Cf-^) (x + m) 



\{^,y) = ACf[Cf-^)c 



9{y>x)- + e{y>x)- 

y y 




P(x) =Cf(Cf-^^ po(l-31nx-ln2x + lnx) +121nx 



+ 



V(x,y) = 



2CF{CF-^)i^Ce{y>x)[F^(^ 
+ ln(l-^)U^ln(l-^) ^ 

^ y^y y ^ 



l_31n^-ln2^ 

y y 



m 




P(x) = CfCa [po - ^ In X - ^ ln2 x) 



V(x,y) = CFCA<Ce{y>x) 



+ 



2y \ yJ\j 



y 2 



(-3 




P(x) =CfCa 



In^ a; + 4 In a; 

+^3o^^ In^ X — Li2(l — x) — Inx Inx^ 



+ 



V(x,y) = CFCAlceiyy x)\^ln'^ f 1--) +- In f 1--) +— Inxln; 

I2y \ yJ y \ yJ yy 



2 y y y 



+ In - In ( 1 



x^ 

y 
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P(x) 



Cf(Cf-^) Po(4Li2(l-x)-ln2x) +81nx 



+ \x 
In |x| ln(l + \x 



V(x,y) = 2Cf(Cf - ^)\ce{y > x)\-\n (l - -) + ^ 
z In - In 1 4Li2 1 +ln^ - 

yy y ^ y^ \ ^ y^ y 



In x In x 

+G'^{x,y) 



q 




P(x) = CfCa Po(i + \t,x + S{x) + ^\t?x- ^Inx + ^In^x) 



-4 In a; - 2hr x 



+ 



V(x,y) = CfCa <j > x) [f'^ ( 2 + 2 In ^ + ^ ln2 ^ + 25 (^) 



In 1 



(l-:^U-ln2fl " 



y 



--In^fl--) --ln(l-- 

y ^ y^ y ^ y 



+ 




P(x) =Cf[Cf 
+CfCa 



Ca\ 



2 J 



Po ( 2 In^ X + 45(x) I - 161nx 



{2 + 25'(x) + ^ In^ X + Inx + ^ In^ x) 



V(x,y) 



Ce{y > x)\f^( 2ln' ^ + 4s(-)] - -In (l - -) 

I- V y ^y^J y ^ y^ 



+CfCa- 



+ 



Ce{y > x)F^ (^4 + 45 (1 - ^) + In^ ^ + 2 In (l - ^) 
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"A 


P(x) = - — CfCa + -C|. + 2CFTrNfjS{l - x) 


V(x,y) = ( - ^CfCa + \cl + 2CFTrN^y{y - x) 


a* b* c* p 

^^^aML^y y/^^^a\y ^^0000^^ y^^^^^ 


P(x) = 

V(x,y) = 



Here S{x) = Li2{x) — Li2(l) and the special notation G^{x) will be clarified in the next section 
(Eq. (3.13)). 

There is, in general, a mixing of quark and antiquark densities in higher-loop calculations. At 
NLO, the diagrams k*, O* contribute to the kernel Piqq expressing the probability to find a quark 
inside a quark (at 1] = 0), they also contribute to the kernel Pigg giving the probability to find 
an antiquark inside a quark. Actually, in the latter case, one should consider two kernels, viz., 
P± — Pigq ± Piqq for 7^ = ±1, [7]. Wc shall separate these contributions and focus on the results 
for Piqg and in the next section. 



3. The structure of the evolution kernels in NLO 

In this section, we discuss the total results of the two-loop calculation and also the general 
structure of the evolution kernels at NLO. We commence with those elements that appear in the 
renormahzation procedure at NLO of both the kernels P and V. 



3.1. DGLAP kernel 



Collecting the "quark-quark" contributions to the NLO DGLAP kernel, presented above, leads 
for Pj^q to the final expression 



where 



Pn(^) 



PLi^) 



Ax - 



Cl ■ Pl{x) + CfCa ■ Pg{^) + CFNfTr ■ Pj,{x) 



Po{x) (3 ln{x) + 4 ln(,x) ln(,T))] ^ + 5(x) (y + 8C(3) - ^) , (3.2) 



(3.1) 



= -2x + 



11 



67 TT^ 



Po(x) \n\x) + — \n(x) + 77 - ^ 



Po{x) \ii{x) + 



9 

+ y S(x) 



,,365 47r2 



4C(3)), (3.3) 
(3.4) 
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This expression together with the expression for Pj^^ can be reduced, after some simple algebraic 
manipulations, to those found in [7]. 

Let us rewrite the expression for P^^^ in Eqs. (3.1)-(3.4) in such a form that corresponds 

to the structure of KiR' at the two-loop level — see Eq. (1.3). Following [12], we consider the 
renormalization of the diagram F and its contracted one-loop subgraph E that can be written 
symbolically as F = ii^ • W , where W is the one-loop remainder. As the result, the pole part of E 
should be multiplied by the finite part of the remainder W and vice versa. All those subgraphs 
Ei that are related to the charge renormalization of the intrinsic vertex in the various diagrams 
(see, e.g. , diagrams d, g, h, 1, m in the list) contribute to the coefficient of the QCD /^-function 



ho 
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^TrNf. After contracting each of these E^ terms, the remainder reduces to one of the 



one-loop diagrams a, b, C from Table 1. The appropriate finite part of each of these diagrams 
in dimensional regularization can be obtained from the differentiation of the auxiliary kernel (cf. 
similar kernels in [12]) P{x;e) — 4x^+^x~^ with respect to the parameter e: 



Po{x) 



e=0 



Po{x) \ii{x) 



(3.5) 



On the other hand, the composite operator illustrated, e.g. . in diagrams a*, b*, C*, O*, q, r 
calls for a different sort of renormalization. Notably, the contracted subgraph E should include 
the composite operator that coincides with that in the one-loop diagrams in Table 1. The latter 
generates the kernel Pq (or Vj^) , while the finite part of the remainder is formed from the finite 
part of -P{x; e), i.e. , by po that finally leads to the contribution 



(po * iPo)+) (x) = po{x) (4 ln(a;) + 4 ln(x) ln(x) — 2 ln^(a;)) . 



(3.6) 



Here the symbol * denotes the Mellin convolution, (/ * g) (x) = dzdy5{x — yz)f(z)g(y). Col- 
lecting all these terms together, one recasts P^ in the form given by the first term in the curly 
brackets below: 



CfPo * 



60I-P0 



+ Po{x)Cf 



67 TT^ 



20 



-NfTr 



+Cf{Cf-^ 
43 



+S{x)Cf 



4x-2{po{x)\n\x)) 
365 _ „ 26 



o.|-c-...r-.(c.-£i)8(c(3)-9 



(3.7a) 

(3.7b) 
(3.7c) 



The second term in the curly brackets in (3.7a) originates from the product of the finite parts of the 
contracted subgraphs Ei, or, more specifically, from the finite part of the charge renormalization 
(diagrams g, h, 1, m) and another finite and specific (see diagrams n*, O*, q, r) part of the 
composite operator, as well as from the pole parts of the remainder that are proportional to Pq. In 
this respect, the coefficient of po appears to be proportional [17] to the two-loop cusp anomalous 
dimension [18], 



^ cusp 



67 



TT 



9 



20 



^ ^ - -^NfT. 
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(3.8) 



The terms in (3.7b) are formed by the diagrams k* and O* with nonplanar elements that also 
contribute to the "quark-antiquark" part P^^^ of the kernel. Finally, the (5-function in (3.7c) 
manifests the fact that the corresponding local current is not conserved. Let us emphasize at this 
point that the expressions in the r.h.s. of (3.7a), (3.7c) together with Eq. (3.8) has the general 
structure of any nonsinglet NLO DGLAP kernel that follows from the renormalization procedure. 
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3.2. ERBL kernel 



Collecting the partial contributions from the diagram-expansion list we arrive at 

V^{x, y) = Cl- V^J(a;, y) + CfCa ■ V^{x, y) + CpN^Tr ■ Vjj{x, y), 



(3.9) 



V^{x,y) = 4C 



X 



X 



e(y > x)- + e{y> x)- + 2G' {x, y) + U C9{y > x) 

y y 



F^ln^- 

y 



+— Inxlnx - -F^ln ^ - (F^ - In - In f 1 - -)] 1 

yy 2 y ^ y ^ y^ ^ j 



+4 



H + 6C(3)-2^ 6{x-y), 



(3.10) 



V^{,x,y) = -2C 



6(y > x)- + 9{y > x)- 

y y 



G^{x,y) 



+ 
+ 



C%>.)2F-(-ln- + --- 
-^-12C(3) + 4^]%-x), 



C^(l/ > x)2F'^(^\n 



X 5 
y 3 



+ 



26,, 
+ y5(y-a;) 



Here and below 



2 

^^(x,?/) = -4c[^(y>x)(F^lnxlny-F^[Li2(x) + Li2(y)] +yF^) +^(y<x) 



X 



(F^-F^)[Li2(l 



(3.11) 
(3.12) 



■In^x] +F^[Li2(j/) -Inxlni/] +F^Li2(x)) . (3.13) 



The term is "diagonal", i.e. , yyG'^{x, y) — xxG'^{y, x), G^ coincides with the similar term G 
in the unpolarized kernel Vi [6] by performing there [10] the replacement F and excluding 

the third term ^F^^ in the first line of (3.13). This term is tied to G^ in order to preserve Fcusp 
in the general structure of . 

The origin of the structure of the ERBL kernel can be considered by analogy with the DGLAP 
explained in [12]. Taking into account that the ERBL auxihary kernel is V{x,y;e) — 
1/2C {9{y > x)P{x/y; e)/y) leads to the following finite part of the corresponding remainder — 
l/2C{9{y > x)po{x/y)/y), which, after introducing an appropriate convolution for the one-loop 
elements [12], (/ g) {x, y) = dzf{x, z)g{z, y), leads to the expression 



V,^{x,y) = J^CFV,''0{bol-V,^)+Ce{y>x)2F^^ + CF[g+, ^^o^j 



X 



+ 5{y - x)Cf 



e{y > x)- + e{y > x)- +2G'{x,y) 
^2 ^ 18 ^ 9 



(3.14a) 
(3.14b) 
(3.14c) 
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The structure of the elements of in (3.14a)-(3.14c) resembles that of in (3.7a)-(3.7c) with 
a natural replacement of notation for the convolution and the symbols (8) — > *, — > po, 
Pq, 2F'^ ^ pq with the exception of the important third term in the curly bracket in (3.14a). In 
addition to the "nondiagonal" term proportional to 60 and to the proper operator renormalization 
(see the first convolution in Eq. (3.14a)), there appears an additional "nondiagonal" term which 
is represented by the commutator [g+, ^Vff] = g+® V(f — V^f ® g+. This term is induced by the 
leading-order anomaly in special conformal transformations of conformal operators, 

9(.,s)=-2C?^^ln(l-i) , (3.16) 

y-x V yJ 

an interesting issue explained in [10, 11]. All the other terms in are "diagonal". Concluding 
these considerations let us mention that the expressions in the r.h.s. of (3.14a), (3.14c) give us 
the elements of the general structure of any NLO ERBL kernel. 



4. Effects of two-loop evolution for the meson DA 



The subject of this section concerns the efi^ects of the two-loop QCD evolution in an appropriate 
example inspired by calculations of the B ^ p ve decay [13, 19]. For the leading-twist DA of the 
transversely polarized p-meson expanded in a Gegenbauer series 



(4.1) 



n=0 



the two-loop evolution of each harmonic from //q to //^, i^nix) — > $„(a;, //^), can be approxi- 
mately represented^ as [12] 



$„(x,/x ) = exp 
Here we have 



''^('^ ) , 7(n,a) 



dr, 



m>n 



7(n, as) = as 7^ {n) + a] 7^(71) , 
(5{as) = - alho- a^hi , 



(4.2) 



(4.3) 
(4.4) 



70 (n) - 70 (m) - 60 



(7o(n)-7o(m)-6o)/6o' 



(n + l)(n + 2) 



4(2n + 3) ""^ 

and 7(n) is the anomalous dimension with 7^ (n) = Z„„/A^„. The coefficients dnm/Nm can be 
calculated analytically (in the form of lengthy sums) by virtue of the knowledge of the structure of 
the "nondiagonal" and "diagonal" terms in expressions (3. 14a)- (3. 14b). Their evaluation for the 
values of the input parameters //q = 1 GeV^, 11% = 36 GeV^ (the latter being the characteristic 
scale of the B ^ p ue decay) for Nf = 4 is presented in Table 2. 



^ The NLO evolution that preserves the renormaUzation-group property for the "nondiagonal" elements was worked 
out in [20]. 
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m\n 





2 


4 


6 


2 


-0.398 











4 


-0.013 


-1.08 








6 


0.024 


-0.297 


-1.269 





8 


0.024 


-0.094 


-0.485 


-1.288 


10 


0.02 


-0.026 


-0.216 


-0.585 


12 


0.015 


-0.001 


-0.103 


-0.303 


14 


0.012 


0.008 


-0.049 


-0.168 


16 


0.01 


0.011 


-0.022 


-0.097 


18 


0.008 


0.012 


-0.007 


-0.057 


20 


0.006 


0.012 


0.000 


-0.033 




Table 2. Numerical values of the dmn/Nm coefficients Fig. 1 The result of one-loop and two-loop evolution 
for the ffi'st 6 harmonics -(An- of (fijl^x; fj,Q — 1 GeV'^) to /i^ = 36 GcV^ 

These coefficients dnm/Nn decrease not so fast as those for the unpolarized case [12]. The numerical 
calculation shows that truncating the sum in (4.2) after the 10th term provides us with a 0.03% 
accuracy at the scale 36 GeV^. 

To get an estimate for the difference between the two-loop and the one- loop result, let us 
compare the corresponding ffist inverse moments of the DA, {x~^) = ip{x)/xdx. The ratios 



) - 1)% at this scale are -4.1% for iio{x), -1.4% for i)2{x), -0.3% for V4(x), 



((^2-loop/ 

and even less for higher harmonics. As regards the model distribution amplitude normalized 

at 



1 GeV^ 



X 



^l) = ^o(x) + 0.29^2(a;) + 0.41^4(a;) - 0.32^6(a;) 



which was obtained obtained for a transversely polarized p-meson in [13], the ratio of the ffist 
inverse moments takes the value 3.6%. The evolution effect on the meson distribution amplitude 
(p'^{x; Hq) is illustrated in Fig. 1. The dashed black line shows the unevolved expression, while the 
result of the two-loop evolution to the scale /x^ is represented by a solid red line, and the one-loop 
result is shown as a blue dashed-dotted hne. 



5. Conclusions 

Let us summarize our findings. In Sec. 2, we presented the diagram-by-diagram results of a di- 
rect two-loop calculation of the DGLAP, P^, and the ERBL, V'^, evolution kernels for transversity 
distributions, employing the Feynman gauge. The mutual correspondence between the V and P 

results, for each of the diagrams, was checked making use of the relation P{z) = lim^^o ^ j , 

[21]. It was found that the total result for P^ coincided with the one in [7] (obtained within a light- 
cone gauge calculation), whereas the total result for turned out to agree with the prediction 
obtained in [10]. 

We worked out the general structure of any nonsinglet NLO DGLAP kernel, Eqs. (3. 7), (3. 8), 
and any NLO ERBL kernel, Eqs. (3.14), respectively, subject to the renormalization procedure. 

The NLO evolution of the DA of twist 2 (for transversely polarized p-meson) was considered 
and its relative effect was estimated for the inverse moment of the corresponding DA. This effect 
amounts to a few per cents (4% for the zero Gegenbauer harmonic) after evolving from the low 
scale /ig — 1 GeV^ to the characteristic scale p^ = 36 GeV^ of the B-decay process. 
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